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A  system  has  been  devised  to  predict  the  neutron- induced  activity 
contribution  to  fallout  exposure  rates.  The  system  uses  the  simplifying 
3  assumptions  of  l)  a  semi- empirical  formula  to  determine  the  soil  capture 

fraction  and  2)  thermal  neutron  cross  sections  to  represent  weapon-neu¬ 
tron  cross  sections.  Results,  using  these  assumptions,  agree  with  those 
of  another  complete  system  for  predicting  the  neutron- induced  activity 
contribution  to  fallout  exposure  rates.  In  audition,  results  from 
various  portions  of  the  system  agree  with  results  obtained  by  more  com¬ 
plicated  methods. 
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SUMMARY 


In  the  general  problem  of  predicting  fallout  exposure  rates  from 
nuclear  detonations  is  embedded  the  problem  of  predicting  exposure 
rates  from  neutron- induced  activity  in  the  fallout.  Certain  restric¬ 
tions,  engendered  by  computer  limitations,  were  necessary  in  order 
that  the  neutron- induced  activity  predictions  be  suitable  for  inclusion 
in  the  solution  to  the  general  fallout  problem.  A  system  based  on 
weapon  neutron  spectra  would  require  a  computer  program  too  extensive 
for  incorporation  in  the  general  program.  Brevity,  without  undue 
sacrifice  of  accuracy,  was  a  requisite.  This  requisite  was  satisfied 
by  assuming  that  weapon  neutrons  could  be  treated  as  thermal  neutrons 
and  that  a  semi -empirical  formula  for  thermal  neutrons  could  be  used 
for  the  purpose  of  making  soil- capture  fraction  and  isotropic- capture 
fraction  predictions. 


The  following  comparisons  each  showed  excellent  agreement  in  the 
region  of  interest:  a)  a  comparison  of  the  thermal  soil-capture  fraction 
predictions  of  this  system  and  a  Monte  Carlo  system  using  weapon-neutron 
spectra;  b)  a  comparison  of  the  predictions  of  the  semi-empirical 
formula  used  in  this  system  with  the  so-called  "exact"  solution;  and  c) 
a  comparison  of  the  exposure  rates  predicted  by  this  system  with  those 
predicted  by  an  extensive  classified  system  using  weapon  neutron  spectra. 
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I.  INTRODUCTION 


T_\ 

A  particle  activity  module-1" '  has  been  developed  at  US1IRDL  to  pre¬ 
dict  exposure  rates  due  to  all  fallout  particles.  The  prediction 
system  discussed  in  this  report  was  devised  to  predict  only  the  neutron- 
induced  activity  contribution  to  these  fallout  exposure  rates.  A  com¬ 
puter  program  has  been  written  for  the  system  by  B.  Hoffman  of  this 
laboratory.  An  adapted  version  of  it  has  been  included  in  the  general 
program. 


Residual  nuclear  radiation  emitted  subsequent  to  an  airburst 
nuclear  explosion  is  almost  entirely  due  to  radioactive  fission  products 
present  in  the  weapon  residues.  However,  if  a  nuclear  weapon  is  deton¬ 
ated  near  a  land  surface,  some  of  the  neutrons  released  in  the  fission 
process  which  escape  the  weapon  will  induce  radioactivity  in  the  soil. 

If  the  point  of  detonation  is  sufficiently  near  the  surface,  the  soil 
will  be  vaporized  with  the  weapon  residues.  Thus,  fallout  particles 
formed  as  the  fire  ball  condenses  will  contain  neutron- induced  activity 
as  well  as  fission-product  activity. 

Although  the  neutron- induced  activity  portion  of  fallout  activity 
is  in  general  quite  small,  there  are  circumstances  under  which  it  may  be 
significant.  It  is  for  this  reason,  as  well  as  for  completeness,  that 
the  neutron- induced  activity  prediction  system  was  incorporated  in  the 
fallout  program.  However,  in  order  to  insure  that  the  neutron  program 
would  not  occupy  a  larger  portion  of  the  generalized  program  than  its 
importance  warranted,  two  major  simplifying  assumptions  were  made: 
l)  all  neutrons  are  emitted  at  thermal  energies;  and  2)  soil  albedo 
for  thermal  neutrons  can  be  determined  semi- empirically.  The  effects  of 
these  simplifications  will  be  amplified  in  the  text;  however,  it  is 
appropriate  to  state  here  that  soil  albedo  (ratio  of  number  of  neutrons 
leaving  soil  to  number  entering  soil)  has  been  found  to  be  nearly  inde¬ 
pendent  of  incident  neutron  energy;  and  also  that,  since  the  soil  is 
volatilized  and  mixed  in  the  fireball,  the  depth  of  neutron  penetration 
before  capture  is  of  no  consequence  (as  it  is  for  prediction  systems 
wherein  the  soil  remains  undist-rubed*). 

This  report  is  arranged  in  such  a  manner  that  the  reader  is  first 
introduced  to  the  rationale  of  the  system.  The  resulting  equations  and 
specifications  of  the  required  input  quantities  follow.  An  estimation 


2)  R) 

See,  for  instance,  Lessier  and  Guy  J  and  Holland  and  Gold''', 


1 


of  the  value  of  each  input  quantity,  based  on  available  unclassified 
information,  is  then  made.  The  associated  computer  program  and  tabu¬ 
lar  inputs  are  relegated  to  the  appendices. 
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II.  METHOD 


As  has  been  noted,  the  neutron- induced  activity  program  was  designed 
to  be  a  part  of  the  generalized  .land  fallout  model  program.  As  such,  the 
form  of  the  output  of  the  neutron  program  was  restricted.  At  the  point 
in  the  general  program  where  the  neutron  program  was  to  be  introduced, 
it  was  required  that  the  neutron  program  output  be  in  the  form: 

total  =  (^i) casing  +  (^l)soil  11X12 °f 
of  roentgens/hr  per  fission/cm2,  where  the  ex¬ 
posure  rates  are  those  found  at  three  feet 
above  an  infinite  plane  uniformly  contaminated 
with  induced  activity.  The  contamination  den¬ 
sity  of  the  induced,  activity  is  expres  »ed  in 
terms  of  fissions  per  square  centimetf  r. 

(Dp) total  is  the  final  output  of  the  neutron  program.  The  general  pro¬ 
gram  adds  this  quantity  to  the  fission  exposure  rates  and  operates  on 
the  sum  to  determine  fallout  exposure  rates. 

The  history  of  weapon  neutrons  from  release  to  capture,  and  the 
conversion  of  the  capture  products  to  induced-activity  exposure  rates 
in  the  required  form  is  traced  in  the  following  section. 

A.  Rationale 

Each  neutron  per  fission  released  in  the  nuclear  process  is  accounted 
for  on  the  basis  of  its  probable  fate.  The  only  possible  fates  are 
assumed  to  be: 

1.  A  neutron  released  by  the  fissioning  of  a  weapon  nucleus 
is  either  captured  by  another  weapon  nucleus,  or  it  is  emitted  by  the 
weapon. 


2.  If  it  is  emitted  by  the  weapon,  it  is  either  captured  by 

the  weapon  casing,  in  which  event  it  is  a  potential  contributor  to  fall¬ 

out  activity,  or  it  escapes  the  weapon. 

3.  If  it  escapes  the  weapon,  it  is  either  captured  in  air  and 
lost*,  or  it  is  captured  in  soil. 

4.  If  it  is  captured  in  soil,  the  soil  is  either  not  volatilized 

and  the  neutron  is  lost,  or  the  soil  is  vola+ilized  and  the  neutron  is 

a  potential  contributor  to  fallout  activity. 

*  The  word  ’".lost’’  here  means  that  the  neutron  is  no  longer  a  possible 
contributor  to  fallout  activity. 
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5.  Those  neutrons  captured  in  the  weapon  casing  or  in  vola¬ 
tilized  soil  are  contributors  to  fallout  activity  if  the  daughter  of 
the  capturing  nucleus  is  radioactive  { i .  e . ,  a  ganma-ray  emitter.  Beta- 
ray  activity  is  not  a  part-  of  this  prediction  system). 

After  finding  the  neutrons  per  fission  captured  in  the  soil  and  in 
the  weapon  casing,  the  fraction  captured  by  each  isotope  of  the  soil  and 
casing  must  be  determined.  Following  this  determination,  the  neutron 
captures  per  fission  by  a  given  isotope  must  be  converted  to  roentgens 
per  hour  per  fission  per  square  centimeter  for  that  isotope.  Finally, 
the  total  exposure  rate  (^v)total  ^oun{*  by  suimning  over  these  isotopic 
exposure  rates. 

B.  Basic  Equations 

The  basic  equation  for  the  induced- activity  exposure  rates  due  to 
activity  induced  in  soil  isotopes  will  first  be  stated  and  then  its 
factors  will  be  examined. 

(Dj)^  (for  soil  isotopes)  =  a  (f  )  JC. 

where  N  ~  neutrons  per  fission  emitted  by  the  weapon 

e~T*~  =  fraction  of  neutrons  emitted  which  escape 

Q .  =  fraction  of  those  entering  soil  which  enter  volatilized 
soil  (solid  angle  fraction) 

a  =  fraction  of  those  emitted  (escaping)  which  are  captured 
in  soil  (soil  capture  fraction) 

(f  1  =  fraction  of  those  captured  in  volatilized  soil  which  are 
1  captured  by  the  i^h  isotope  of  the  soil  (isotopic  capture 
fraction) 

K.  =  roentgens  per  hour  due  to  }S0tope  per  it**1  isotoee  atom  de¬ 
posited  per  cm2.  The  exposure  rate  bcirg  for  a  point  at  3 
ft  above  an  infinite  plane  on  which  one  nucleus  of  the  jth 
type  is  deposited  per  cm2,  h'ote:  K.  is  time  dependent. 

The  first  five  factors  simply  determine  how  many  neutrons  per 
fission  are  captured  by  the  i^h  isotope.  Ihe  final  factor,  Kj,  then 
converts  this  to  exposure  rate  in  the  appropriate  form.  The  equivalent 
equation  for  casing  isotones  is 

(D_).  (for  casing  isotores)  =  N  (X-e~  ^)  (f  ).K. 

\  JL/1  v  O  C  l  Z 

where  (l-e  )  =  casing  capture  fraction 

-  When  the  casing  is  assumed  to  be  r*f  zero  th/ehness, 
is  equal  to  the  lumber  cscaninf.  fee  Section  C-V. 


the  number  emitted 


The  basic  equation  for  the  total  neutron- induced  activity  contribu¬ 
tion  to  fallout  exposure  rates  is  then 

(D  )total  -  7.  (D_) .  (summed  over  soil  isotopes)  + 

I  i  1  1 

+  E  (D^)^  (summed  over  casing  isotopes) 

i  x 

C.  Estimation  of  Input  Values 

1.  Neutrons  Emitted  oer  Fission,  T?  . 

4  \  o 

(see  E&U  ,  Sections  1.4l  -  1 ,6b)* ** 

a)  From  1.43,  " ...when  a  nucleus  captures  a  neutron  and  suf¬ 
fers  fission  f  neutrons  are  released;  let  l  be  the 
average  number  of  neutrons  lost .... 

b)  From  1.44,  "For  uranium-235,  f  is  about  2.5,  t  may  oe 
taken  to  be  roughly  0.5. -  - 


Since  it  takes  one  neutron  to  cause  a  nucleus  tc  fission,  there  is 
an  overall  excess  of  2.5  -  1  =  1.5  neutrons  per  fission  emitted  by  the 
fission  process.  Therefore,  for  a  pure  fission  weapon. 


=  1.5 


fission  neutrons 
fission 


Note  that  this  is  an  overall  excess  and 


includes  all  of  those  neutrons  released  in  the  last  generation  of  fission 


as  well  as  the  i  =  0.>  neutrons  per  fission  lest  during  the  chain  reaction. 


For  a  fission-fusion  weapon,  we  use  the  follcwing  to  estimate  the 
number  of  fusion  neutrons  emitted  per  fission: 

c)  from  1.66,  "H2  +  II2  =  Ke3  +  n  +  3.2  Mev 

-22^1, 

h  +  H  =  H-'  +  II  +  4.0  Mev 

H~  +  H2  =  HeU  +  n  +  17.6  Hev-..." 

or  pH*  -H1  +  He3  +  +  2n  +  24.8  Mev. 

This  shows  that  two  neutrons  are  released  per  24.8  !*ev 
(fusion). 


*  The  practice  of  using  EKV?,  or  any  oilier  source  of  reliable  unclassified 
information,  is  maintained  throughout  this  report  to  deduce  reasonable 
input  values.  Naturally,  the  classified  weapon  data  for  a  specific 
weapon  are  the  correct  input  values  to  be  used. 

**  Note  that  "lost"  here  means  lost  to  the  fission  or  fusion  process. 


d)  Since  2.625  x  10^  Mev  is  equivalent  to  1  kiloton, 

2  neutrons  2.625  x  10^  Mev  0  ^  24  fusion  neutrons 

54.6  Mev  x  KT  ikilcton)  =  *  x  1U  KT  (fusion) 

e)  Pro©  1.4l,  1,45  x  10  J  fissions  is  equivalent  to  1  KT 
(fission);  and 


d  t/f  = 


0.5  _  1  neutrons  lost 
2,5  "  5  neutrons  released 


On  combining  c),  d),  e),  and  f),  ve  find 


„  /fusion  neutron  s\  1 
o\  fission  /  5 


neutrons  lost 
neutrons  released 


0-10  ,^24  (fusion  neutrons  released] 

C.lc  X  1U  - ^  T„ . ...  ■  .. ,  « 

KT  (fusion) 


1.45  x  10 


fissions 
KT  (fission) 


VI  [KT( fusion)]  ^  , 

_F _ „  r  fusion  neutrons 

u  _ •  \t  -  **  fission 

1 KT( fission f 


The  total  number  of  neutrons  emitted  per  fission  by  the  weapon  is  then 

^  ,  -  fission  neutrons  +  fusion  neutrons 

*0  fission  fission 

No  prevision  has  been  made  to  calculate  NQ  from  this  formula  in 
the  computer  prog  ram.  1?^  itself  is  the  input  quantity. 

2.  Casing  Capture  Fraction,  (l-e~^). 

In  the  version  of  this  prediction  system  incorporated  in  the  gener¬ 
alized  la;>d  fallout  model  computer  program,  the  casing  capture  fraction 
has  been  omitted.  The  quantity  usually  reported  in  the  classified  lit¬ 
erature  is  the  number  of  neutrons  escaping  the  casing.  If  this  quantity 
is  used  for  5  ,  the  casing  capture  fraction  is  automatically  zero.  Thus, 
since  compute?  space  is  at  a  premium  in  the  general  program,  the  reten¬ 
tion  of  the  ability  to  account  for  the  casing  effects  may  be  unjustified. 
However,  e  method  for  calculating  casing  captures  has  been  retained  in 
the  program  reported  here,  since  computer  space  is  not  at  a  premium  and 
scae  estimation  of  the  casing  effects  cay  prove  desirable. 


*  That  is,  ve  assume  the  ratio  of  neutrons  "lest’*  to  neutrons  released 
for  the  fission- fusion  case  to  be  essentially  that  given  by  Section 
1.44  of  n«W.  Higher  or  lower  fractions  could  be  inserted  at  this 
point  as  desired  by  a  user. 
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In  the  expression  ( * — c  ),  £  is  the  integral  macroscopic  crosr- 

section  and  X  is  the  effective  casing  thickness,  i.e., 

l5A  .  -i 

£  -  o—  o  i  i  ntegra-1 )  cr 

where  n  =  casing  density 

Avogadro's  nun’.  r 

A  =  effective  atomic  weight  of  casing  elements 
and  7  interval  capture  cross  section 

Effective  casing  thickness  means  the  oath  length  the  neutrons 
must  travel  in  casing  materials.  The  determination  of  this  length  re¬ 
quires  a  thorough  knowledge  of  all  casing  parameters  and  of  the  angular 
and  energy  distribution  of  neutrons  incident  on  the  internal  face  of  the 
casing.  The  determination  of  X  must  be  made  independently  of  this  pre¬ 
diction  system.  The  same  is  true  of  a  (integral),  for,  as  noted  in  the 
introduction,  all  neutrons  are  emitted  at  thermal  energies  in  this  pre¬ 
diction  system.  Thus,  the  computer  calculates 

,  *A  »A 

£  =  p  -r-  a  (thermal),  not  £  =  p  ~  o  (integral) 

A 


To  correct  for  this,  an  adjusted  effective  casing  thickness, 

/  y  °  (integral) 

‘  -o  (thermal) 

is  entered  as  input  so  that 
z'x'  -  £X 


and  the  desired  numerical  value  is  achieved.  The  computer  requires  as 
input  the  casing  density  and  fraction  by  weight  of  casing  elements  to 
calculate  £* .  The  neutrons  per  fission  captured  by  the  casing  are 
then  treated  in  the  same  manner  as  those  captured  in  soil,  as  will  be 
discussed,  to  arrive  at  the  exposure  rates  from  casing- induced  nuclides. 


To  simplify  the  discussion  in  the  balance  of  this  report,  no 
distinction  hI  J  <  be  ’-a-ic  between  the  ton's  "icu:  ronr.  emitted  and 

escaping ’ - 


f>trther 

p.cutroi 


o* 


Solid  Angle  Fraction, 


In  Figure  1,  h  <  <  one  capture  mean  free  path  for  neutrons.  From 
the  figure,  it  is  seen  that  all  neutrons  emitted  from  the  weapon  in  a 
downward  direction  enter  the  soil.  This  is  equivalent  to  stating  that 
the  earth  subtends  a  solid  angle  of  2 it  steradians  and  t  hat  air  attenua¬ 
tion  is  zero.  For  the  heights  of  burst  of  significance  both  statements 
are  nearly  exact. 


Neutrons  entering  this  region  ore  lost 
=  ’4  of  those  emitted 


BURST  POINT 


Fig.  1  BURST  GEOMETRY 


If  we  denote  the  solid  angle  subtended  by  the  volatilization  crater 
as  Q,  and  the  ratio  of  the  solid  angle  subtended  by  the  crater  to  that 
subtended  by  the  earth  as  Q^,  then  this  ratio  (solid  angle  fraction)  is 

Since,  in  polar  coordinates,  the  crater  solid  angle  can  be  expressed 


Q  =  2r,  J  sin  9d  9  =  2n  (l-cos  3  ) 

o  max 


8 


and,  as  seen  from  the  figure. 


cos 


G 

max 


h 


we  find. 


In  order  to  determine  fit  for  a  given  height  of  burst,  h  (ft.),  and 
weapon  yield,  W  (KT),  we  must  determine  the  crater  radius,  Cjj  (ft.),  as 
a  function  of  h  and  W.  In  order  to  find  this  function  we  again  turn  to 
ENlf4)  as  a  convenient  source  of  unclassified  information: 

a)  from  6.08,  "It  has  been  estimated  that  for  a  1-kiloton 
nuclear  burst  near  the  surface,  the  diameter  of  the  crater, 
i.e.,  of  the  hole,  will  be  about  130  feet  in  dry  soil...." 

b)  From  o.Oy,  "...  for  an  explosion  of  U  kilotons  yield.,  -  on 
the  surface,  the  diameter  ...  of  the  crater  will  be  If  times 
the  values  quoted  above  for  a  1  kiloton  burst." 

c)  From  2.117,  ("  .  ..R  is  the  fireball  radius  in  feet  and  W 
is  the  explosion  yield  in  kilotons...."') 

o  U 

"R  (at  breakaway)  for  air  burst  ~  110  V  *  " 

"The  size  of  the  fireball  is  not  well  defined  in  its  later 
stages,  but  as  a  rough  approximation  the  maximum  radius  may 
be  taken  to  be  about  twice  that  at  t«c  time  of  breakaway. . . 

d)  From  6.11,  "As  the  height  of  burst  increases,  the  dimensions 
of  the  center  vary  in  a  rather  complicated  manner. ...  In 
fact,  for  an  appreciable  crater  to  be  formed,  the  height  of 
burst  should  be  not  more  than  about  one- tenth  of  the  maximum 
fireball  radius.” 

e)  From  2.05,  "...the  fireball  from  a  1-megaton  weapon  ...  in¬ 
creases  to  a  maximum  value  of  7200  feet  (across)." 

Using  a  and  b,  we  find 

.1 

C  ~  65  W  at  h  =  0. 

R 


Q 


Now,  if  we  use  (from  c  above)  R  (at  breakaway)  for  air  burst  «  llOtf^*^ 

and  R  (maximum)  ~  2R  (at.  breakaway)  and  combine  these  vit-h  C  =0  for 

h  ^  V;  maximum  fireball  radius  (from  d  above),  we  find  C_  =R0  at 

K 

h«r  V-  x  2R  (at  breakaway) 

«  M  x  2  x  now0*11 

»  22N0*** 

0  U 

Solving  h  ss  22W  *  for  a  1-megaton  burst,  we  find  CR  =  0  for 
h  ~  349  ft.  and  W  =  1  2-tF.  Alternatively,  if  we  use  (from  e  above) 

R  (maximum)  =  =  36 00  ft  for  W  1  MT, 

then  CL«0athas3^x  36  00  =  3^0  ft.  We  wish  now  to  test  the  effect 
of  assuming  cube  root  scaling  for  height  of  burst.  For  a  1-megaton  burst, 

1 

«  0  at  h  =  36w5=  360  ft. 

x  lx 

Thus,  h  =  3&W  is  within  a  few  percent  of  h  =  22W^*  at  1-KT.  It 

agrees  exactly  at  about  1.6  MT  and  remains  within  1*0 £  from  1  KT  to  100 
r?T.  Since  all  of  the  above  information  from  EEW  is  approximate,  it 
seems  reasonable  to  use  cube  root  scaling  fcr  height  of  burst,  since  its 
use  simplifies  the  calculations. 


Since  the  crater  radius  decreases  in  a  rather  complicated  (and 
unknown)  manner  as  the  height  of  burst  increases,  the  simplest  expression, 
i.e.,  linear,  was  used  to  relate  crater  radius  to  height  of  burst.  Thus, 

n  .  65Hi  .  Si!  h 

R  22W3 

or,  letting  CR  =  CRA3  and  h  =  h/W3, 


S' _  ^ 

C  =  65  -  1.6  h 

Inserting  this  result  into  the  expression  for  Q^.  we  have 
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(h  >36) 


The  computer  input  quantities  are  h  (ft)  and  M  (KT).  The  general  program 
tests  for  h  >  36;  however,  this  program  does  not.  The  program  should  not 
he  used  for  h  >  36. 

4.  Soil  Capture  Fraction,  a 

The  neutron  number  albedo  or  reflection  coefficient  of  soil  may  be 
defined  as 

A  =  neutrons  exiting  soil/neutrons  entering  soil. 

For  an  infinite  medium,  the  neutrons  are  either  captured  or  reflected 
therefore, 

1  -  A  =  neutrons  captured  in  soil/neutrons  entering  soil. 

The  soil  capture  fraction,  a,  is  defined  as, 

a  =  neutrons  captured  in  soil/neutrons  emitted  by  the  weapon 

Using  the  geometry  assumed  in  Fig.  1,  it  is  seen  that 

\  =  neutrons  entering  soil/neutron  emitted  by  the  weapon; 

therefore, 

c  =  i(i-A) 

Kalpem^  has  derived  a  theoretical  formula  for  the  neutron  number 
albedo  of  the  plane  surface  of  an  infinitely  thick  medium,  given  an 
incident  isotropic  distribution  of  thermal  neutrons.  The  formula  may 
be  written 


where  o  *  and  a  *  are  the  respective  thermal  neutron  scattering  and 
SSC  s 

capture  cross  sections  of  a  given  soil.  Using  iialpern’s  formula, 
a  (thermal)  =  jt(i-A) 


In  order  to  determine  the  effect  of  the  simplifying  assumptions  on  ,, 
the  soil  capture  fraction  predictions,  the  Konte  Carlo  values  of  Biggers 
were  used  to  calculate  soil  capture  fractions  for  each  of  48  energy- range 
groups.  The  results  are  shown  in  Table  1. 

*  The  methods  of  calculating  a  and  c  are  given  in  Section  C.5- 

ssc  s 


TABLE  1 


Range 

\(  meters) 
Energy  n. 

(«ev) 

SOIL  CAPTUR 

(Biggers 1 

0-250  250-500 

FRACTIONS 

Soil) 

5  CO- 1000 

>1000 

a  (E) 

.0005 

.118 

(.50) 

.139 

(.094) 

.159 

(.002) 

(.000) 

.120 

(.010) 

.001 

.121 

(.90) 

.132 

(.098) 

.155 

(.002) 

(.000) 

.122 

(.107) 

.01 

.116 

(.82) 

.126 

(.170) 

.107 

(.010) 

(.000) 

.117 

(.107) 

0.1 

.111 

(.73) 

.116 

(.240) 

.122 

(.030) 

(.000) 

.112 

(.107) 

1 

.100 

(.58) 

.111 

(.300) 

.107 

(.070) 

.077 

(,o4o) 

.103 

(.280) 

3 

.151 

(.53) 

.118 

(.310) 

.102 

(.160) 

.106 

(.on) 

.126 

(.1263) 

4 

.138 

(.51) 

.115 

(.310) 

.116 

(.160) 

.108 

(.020) 

.127 

(.1263) 

6 

.172 

(.51) 

.129 

(.310) 

.130 

(.150) 

.123 

(.020) 

.151 

(.1264) 

8 

.193 

(.53) 

.l4o 

(.280) 

.128 

(.160) 

.123 

(.030) 

.166 

(.0025) 

10 

.191 

(.54) 

.148 

(.270) 

.131 

(.160) 

.117 

(.030) 

.168 

(.0025) 

12 

.188 

(.54) 

.l4l 

(.270) 

.124 

(.170) 

.104 

(.030) 

.163 

(.0025) 

l4 

.188 

(.53) 

.130 

(.270) 

.125 

(.150) 

.126 

(.030) 

.160 

(.0025) 
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The  entry  without  parentheses  in  each  bin  is  the  soil  capture 
fraction  a( R,E)  for  Riggers*  soil  for  neutrons  emittted  at  energy  E  and 
entering  the  soil  in  the  range  R.  The  entry  in  parentheses  is  the  ratio 
of  the  number  of  neutrons  entering  the  soil  in  the  given  energy-range 
bin  to  the  total  number  emitted  at  the  given  energy.  These  parenthetical 
entries  are  used  as  range  weighting  factors,  vR,  to  determine  each  ct(E) 
shown  in  the  final  column  of  the  table.  The  entry  with  parentheses  in 
the  c*(e)  column  is  the  ratio  of  the  number  of  neutrons  of  the  given 
energy  to  the  total  number  of  neutrons  emitted  by  a  fission  weapon  (a 
Watt  fission  spectrum  was  usea)^'.  These  parenthetical  entries  are  used 
as  energy  weighting  factors,  wR,  to  determine  a  (fission  spectrum). 

Thus, 


and 


a(E)  =  L  wR  a(R,E) 
R 


a  (fission  spectrum)  =  £  wR  a(E). 

E 


•Hie  result  is 


a  (fission  spectrum)  =  .120  neutrons  captured/neutrcn  emitted. 

When  a  (thermal)  is  calculated  using  the  formula  based  on  Halpem, 
with  the  a3  and  oSSc  being  calculated  for  the  soil  used  by  Biggers,  the 
result  is 


a  (thermal)  =  .119 

Thus  the  a  ( thermal)  formula  agrees  with  the  calculations  based  on 
Biggers*  data.  Further,  an  inspection  of  the  a(E)  column  shows  the 
capture  fraction  to  be  highly  insensitive  to  the  energy  of  the  emitted 
neutron  so  that  spectral  differences  should  not  affect  the  results. 

Another  investigator,  Rafalski^,  reports  the  solution  of  the  so 
called  "exact"  equation  for  the  reflection  of  a  parallel  beam  of  neu¬ 
trons  normally  incident  on  an  infinitely  thick  wall. 

The  equation  is 

*  -  *  -  -r  a  -  *)  Vi  -  %  '(rfL, 

where  £^  and  T.r  are  the  macroscopic  total  and  scattering  cross  sections, 
the  function  is  found  from  the  integral  equation 
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and  the  coefficient  k  is  found  from  the  transcendental  equation 
tanh  k  —  =  1- 


In  Appendix  A,  the  sample  output  shows  a  soil  for  which  crg  =  .203* 
a  =9.16  and  the  resulting  a  =  .170.  When  these  cross  sections  are 
uli&  in  conjunction  with  the  "exact"  solution,  a  =  .Xj6  (graphical  inter¬ 
polation  was  used  and  the  last  figure  is  uncertain) .  In  still  another 
comparison  {using  the  cross  sections  Biggers  reported  for  his  soil  ) 
d  -  .145  (" exact"- interpolated)  while,  using  the  same  cross  sections 
with  the  prediction  system  formula,  a  =  . l4o. 

All  intercomparisons  in  this  section  indicate  that  the  simplifying 
assumptions  are  quite  justified  for  soil,  and  that  the  simple  soil  cap¬ 
ture  fraction  formula  used  in  this  prediction  system  is  acceptable. 

The  computer  has  been  programmed  tc  compute  a,  a  for  a 

given  soil  given  the  fraction  by  weight  of  each  element  in  the  soil,  the 
abundance  of  each  isotope  of  the  element,  and  the  scattering  and  capture 
cross  sections  of  each  isotope.  The  details  are  given  in  the  following 
section. 


5.  Isotopic  Capture  Fraction,  (f  ). 

Given  that  a  neutron  is  captured  in  the  soil,  we  wish  to  knew  the 
probability  that,  it  is  captured  by  a  given  isotope  of  a  given  element 
of  that  soil. 

Let  f^  =  number  of  atoms  of  i—  isotope  per  atom  of  soil 

=  thermal  neutron  capture  cross  section  of  i—  isotope 


Then,  since  q-  =  thermal  neutron  capture  cross  section  of  soil  (as 
previously  defined 


and,  since  (f  ).  =  neutrons  captured  by  i—  isotope  per  neutron  captured 
in  soil  (as  previously  defined). 


Ik 


t  +  \ 

v*c'i 


[Neutron  resonances  in  light  nuclei  (Z  <  30}  are  almost  completely 
scattering0',  and,  since  soil  elements  are  nearly  all  in  this  category, 
they  are  presumed  to  follow  closely  the  l/v  law;  therefore,  the  relative 
captures  are  proportional  to  thermal  cross  sections  and  no  integral 
cross  section  adjustment  for  weapon- spectrum  energies  need  be  used  to 
refine  the  isotopic  capture  fraction  estimate.] 


The  values  of  f.  are  not  generally  available  in  the  literature.  The 
fraction  by  weight  ol  the  kit  element  in  a  given  soil  (f  1,  is  the 
quantity  usually  reported.  As  a  consequence,  provision  nab  been  made  in 
the  computer  program  to  convert  this  quantity  to  f^  by  the  relationship 


fi  = 


(fvV\  ,  . 


where  A.  =  atomic  weight  of  k —  element 

&  till 

and  f.,  =  atoms  of  i —  isotope  per  atom  of  k. 


..  The  program  also  obtains  A. ,  given  f..  and  A 
isotope  of  kit  element  by  K  ix.  : 


atonic  weight  of 


*k  *  f  f  ik  Aik 

The  computer  inputs  are  (f^)^,  f^  and  at.  A  tabulation  of  the 

(f  )  values  of  18  elements  in  three  soils  is  contained  in  Appendix  B, 
TaSle  1.  The  f.  and  a  values  of  the  53  isotopes  of  these  elements 
are  contained  in  Appendix  3,  Table  2.  Also  contained  in  Table  2  are 
the  scattering  cross  sections  of  the  isotopes,  (®sc)j»  used  as  input 
to  calculate 


ssc 


=  Z 
i 


ri  (V>i 


ibis  is  the  soil  scattering  cross  section  used  in  the  program  to  cal¬ 
culate  the  soil  capture  fraction,  a.  The  cross  sections  and  isotopic 
fractions  are  taken  from  BNL  3^5°)  and  from  the  G.  E.  Chart  of  the 

Nuclides^). 


^yjUpjr-j 
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6.  Conversion  to  Exposure  Rate,  k 
Referring  back  to  the  basic  equation 

<Vi  ■  v'“  V  <f A  K1  - 

the  regaining  quantity  to  be  discussed  is 

[B/hr  j  neutron  captured  by  i—  isotope/cn2] . 

th 

For  each  neutron  captured  by  a  nucleus  of  the  i —  isotope,  a 
daughter  nucleus  is  produced.  He  are  concerned  here  with  radioactive 
daughters,  each  of  which  has  a  characteristic  decay  scheme  with  its 
own  decay  constant. 


and  its  own  set  cf  characteristic  gaana-rays  and  their  emission  ratios, 
(fy)i1  [photons  of  energy/ disintegration  of  stow  of  iH  isotope] . 

Let  R 

1 

the  conversion  factor  relating  the  gaae*-ray  exposure  r^ge  at 
3*  above  an  infinite  plane  to  the  emission  of  one  photon  of  the  J — 
energy  per  second  from  each  square  centimeter  of  that  plane.  Then 

Ki  *  V  v  *  (Vi jRj 

where  t  is  the  tine  after  capture  (Tine  after  capture  is  taken  to  be  tine 
after  detonation). 


[R/hr  j  photon  of  energy  per  as  per  sec]  = 


\ppendix  B,  Table  2  contains  one  additional  quantity,  (f  1 ).,  the 
half-life  of  the  daughter  of  the  iHi  isotope.  Shis  is  the  quintity 
usually  reported  in  the  literature.  Provision  bos  been  aade  in  the  com¬ 
puter  program  to  convert  it  to  X..  by 


.  -693 


and  to  convert  each  Ti  ,  in  whatever  standard  tine  units  it  has  been 
reported  (seconds,  niftutes,  hours,  days,  years),  to  X  in  seconds. 
Therefore  t  is  entered  in  seconds. 
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Th'.-  decay  schemes  of  the  isotopes,  -ilso  taken  frca  Reviews  of 
Modern  Physics,10'  yielded  the  (f^). .  for  each  isotope.  They  are  listed 
la  Append i*  3,  Table  3*  The  only  nor.-sta.  lard  quantity  used  in  the  con¬ 
version  was  R  «.  The  calculations  were  made  by  S.  C.  Rainey*  of  this 
laboratory  using  the  Cates  end  Eisenhauer*-1'  build-up  factors.  They 
are  also  tabulated  in  Appendix  B,  Table  3. 


V 


The  inputs  to  the  computer  for  the  conversion  are  (T* ).,  t,  and 

2  1 


♦Private  ccoB&unicatioo 
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HI.  DISCUSSION  AMD  CONCLUSIONS 


We  wish  to  distinguish  between  the  internal  accuracy*  of  the  predic¬ 
tion  system  and  the  overall  accuracy.  Internal  accuracy  means  the  accuracy 
of  the  assumptions  used  in  calculating  a,  v (f  )  ,  and  Q_.  H  ,  EX,  and  K. 
are  external  calculations,  i.e.,  the  data  baSe1^  input  to  °the  system  1 
and  would  be  required  for  any  other  system. 


The  only  internal  part  of  (fc).  is  the  use  of  cr.  (thermal)  rather 
than  the  unknown  o.  (integral),  which  is  variable,  not  only  with  weapon 
neutron  spectra,  bat  with  depth  in  the  soil.  The  evidence  advanced  to 
verity  the  formula  fora  seems  to  indicate  that  the  use  of  thermal  cross 
sections  is  Justified  for  the  soil  capture  fraction;  therefore,  any 
biasing  in  the  soil  capture  fractions  due  to  the  use  of  thermal  orosa 
sections  would  appear  to  be  considerably  less  than  any  biasing  due  to 
erroneous  (fv)^  values. 


It  may  be  noted  that  the  high-crosc  section,  trace  elements  Gd,  Sm, 
end  B  have  been  included  in  the  (fV  tabulations  of  Appendix  B,  Table  2. 
To  show  the  effect  of  their  presence?  we  let 


neutrons  oroducd 


ng  radioactive  nuclides 
fission 


then 

p  .(oafa£  (fc)i 
active 


where 

£(f  ).  *  the  isotopic  capture  fractions  sussed  over  radioactive 
active  daug?lters 
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we  can  write 


S  = 


1.155  qji  2;  r.  o. 

active 

+ ° r?  "<o* 


ssc 


S' 


The  neutron- induced  daughters  of  Gd,  Sm,  and  B  isotopes  are  stable; 
therefore,  the  summation  over  the  active  isotopes  is  affected  by  the 
presence  of  these  trace  elements  only  to  the  extent  that  the  isotopic 
fraction  of  each  active- daughter  isotope  is  decreased  by  the  inclusion 
of  Gd,  So,  and  B  in  the  soil.  Since  the  elements  appear  in  soil  in 
minute  concentrations,  the  numerator  is  insensitive  to  their  presence. 
The  same  argument  holds  true  for  the  <jssc  in  the  denominator  since 
these  trace  elements  have  negligible  scattering  cross  sections.  However, 
the  capture  cross  section  o&  does  Increase  somewhat  by  their  inclusion, 
and  9  will  go  down  roughly  as  (crs)~2  increases  (in  general  o  »o  ). 

To  give  a  quantitative  example  ve  use  c^=  .203  from  the  sample  out¬ 
put  in  Appendix  A.  This  a  was  calculated  using  the  STS  (Folan)  soil 
listed  in  Appendix  B,  Table  1.  For  this  soil  the  (f  )  of  Gd  =  0.  Let 
us  assume  the  true  Gd  content  should  have  been  that  listed  for  the 
Earth's  Crust,  Appendix  B,  Table  1,  (f  L  *  6.4  x  10"°.  To  simplify 
the  calculation  we  use  the  atomic  weight  and  cross  section  of  the  ele¬ 
ment,  as  taken  from  the  Chart  of  the  Kuclides?/,  rather  than  make  the 
calculations  based  upon  the  individual  isotopes  of  the  element.  From 
the  chart  then. 


Atomic  Weight  of  Gd,  A^  =  157.25 


and 


=  4.6  x  10  bams/atoms  (thermal  cross  sections) 


~-2 


A  hand  calculation  shows  E(fw)k/A^  =  7.2  x  10”  for  HIS  (Folan). 


Therefore,  using  A^  to  get  the  atomic  fraction,  f^, 

_  ^jpGd/^Gd 
fcd  -  z 

_  6.4  x  10~6/157.25 
7.2  x  10~2 


=  5.65  x  10"7 


_ atoms  of  Gd  _ 

atom  of  NTS  (Poian)  soil 
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low,  remembering  that 
"1 

*  .203  (without  Gd) 
When  Gd  is  added,  ve  hare 


-203  +  fGd°Gd 

.203  ♦  5*65  x  IO"7  x  4.6  x  10* 
.203  +  *026 


a  .229  hams  (with  Gd), 


so  that  the  addition  of  6.4  parts  per  Million  (by  weight)  of  Gd  increases 
e  hy  aore  than  lOjt.  Bie  ratio  of  the  p's  then  is 

g(vith  Gd)  /  .203 

pfvithout  Gd)  "y.229 

»  .942 

and  the  change  is  Minor. 


If,  however,  a  soil  analysis  is  only  Made  for  elements  having  an 
abundance  of  at  least  10 T*  parts  hy  weight,  then  as  Much  as  10"”4  parts 
of  Gd  could  he  present  without  being  reported,  then  f_.  *  8.83  x  10”° 
and  0,1 


o  *  .203  ♦  8.83  x  10"6  x  4.6  x  10* 
8 

*  .203  ♦  .4o6 
a  .609  (with  Gd) 

and  the  ratio  of  the  P's  becomes 


*  .58 

Thus,  neglecting  10~*  parts  of  Gd  would  result  in  a  significant  over- 
protection  of  exposure  rates. 
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If,  now,  ve  investigate  changes  in  0  as  a  function  of  changes  in  the 
(fv)  values  of  the  significant  contributors  to  exposure  rates  of  the 
sample  output  in  Appendix  A,  ve  see  that  Na(Z*ll)  and  Mn(Z=25)  contribute 
better  than  of  the  exposure  rate;  therefore,  the  change  in  3  is 
nearly  directly  proportional  to  any  change  in  the  weight  fractions  of 
these  elements . 

The  weight  fractions  of  the  significant  contributors,  like  those  of 
the  trace  elements,  are  of  course  external  to  the  system,  and  like  Hq,  LX, 
and  would  be  required  by  any  csher  system,  regardless  of  its  sophisti¬ 
cation. 

With  respect  to  the  solid  argle  fraction  should  be  correct  for  a 
surface  burst  since  no  calculations  using  h,  Cjg,  and  V  are  involved.  The 
greatest  error  in  will  be  for  low  yield  weapons  detonated  near  the 
maximum  height  of  burst.  But,  in  this  region,  the  per-cent  error  becomes 
academic  because  the  neutron-induced  activity  contribution  to  fallout 
approaches  zero  as  height  of  burst  Increases  while  the  fission  contribu¬ 
tion  remains  constant. 

TO  check  the  internal  accuracy  of  our  system,  its  output  for  given 
conditions  was  checked  against  a  similar,  but  more  elaborate,  system  re¬ 
ported  in  the  classified  literature  by  Polan  et  al.  Polan *s  input  was 
the  spectral  distribution  of  neutron  energies  for  each  of  six  weapons 
with  varying  fission- fusion  ratios.  The  total  number  of  neutrons  reported 
as  being  emitted  by  each  weapon  was  converted  to  H  for  use  in  our  pre¬ 
diction-system.  All  weapons  were  caseless  and  detonated  at  surface  zero, 
thus  e  and  Q_  were  not  Involved.  The  average  of  the  ratios  of  the 
predictions  for  the  six  weapons  showed  Bolan's  to  be  0.8£  higher,  with 
the  maximal  disagreement  showing  Polan  to  be  2.7 f»  higher. 

Some  of  the  elements  Included  in  our  system  were  not  included  in 
Folan's  and  vice  versa;  therefore,  to  Insure  that  the  agreement  was  not 
wholly  fortuitous,  a  hand  calculation  was  made  of  o  and  a  using 
Polan* s  soil  composition.  From  these,  the  soil  capture  fraction,  a, 
was  calculated  and  compared  with  Polan* s  only  reported  capture  fraction 
for  one  soil  and  one  weapon.  The  ratio  of  the  capture  fractions  showed 
Polan *8  to  be  k.k%  higher. 
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Modules  for  a  Generalized 

2 

and 


Fail 


Doses  fron 


fcutron-Induced  Besldual  Radioactivi 


11.  L.  D.  Gates,  Jr.,  C.  Eisenhauer,  Spectral  Distribution  of  Gamma 
Rays  Propagated  in  Air,  Armed  Forces  Special  Weapons  Project, 
Technical  Analysis  Report,  AFSWP  Ho.  502A  (195*0  • 

12.  M.  Polac,  et  al..  Contributions  to  Fallout  from  Neutron  Activated 
Soil  (u).  Ford  Instrument  Company,  DASA  1562,  December  1964, 

"  (Classified)" . 
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APPENDIX  A 

Computer  Program 

Equivalent  Symbols  and  Gloaaaxy 

Program 

Sample  Output 
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Equivalent  Symbols  of  Text  and  Program 


Text 

PmffruiK 

Kc>i 

SIGISC 

IW(L) 

fik 

fai(i,l) 

°i 

sigi(i,l) 

T* 

HL 

'Vij 

FOG(N,I,L) 

RJ 

RNY(N,  I,  L) 

°s 

SIGS 

Di 

SDBL 

H 

WYK 

h 

HOB 

X 

THICK 

p 

HHOC 

Bo 

EKEBf 

*** 

h 

SCBOB 

°f 

FOM 

ac 

S3GC 

E 

SIGH 

a 

ALB 

o 

8SC 

SIGSSC 
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FORTRAN 


31 79  SYSTEM 


IHOOCEO  (MAIN) 


FORTRAN* 


C  12  JAN  66  3179  SYSTEM  -  INOOCEO  JONES  /  HOFFMAN 


C  GLOSSARY  - 

C  ALB  -  NEUTRONS  CAPTURED  RV  SOIL  PER  NEUTRON  EMITTED 
C  EMITN  •  NEUTRONS  EM1TTEO  PER  FISSION  OF  THE  WEAPON 
C  FAllItL)  •  ATOMS  OF  ISOTOPE  I  OF  L  PER  ATOM  OF  L 
C  FAT(L)  »  ATOMS  OF  L  PER  ATOM  OF  SOIL 

C  FOGIN# I  *L)  *  PHOTONS  OF  A  PARTICULAR  ENERGY  PER  DISINTEGRATION 
C  FOM  •  SOLID  ANGLE  FRACTION  SUBTENDED  BY  VOLALITI2ATION  CRATER  OF 
C  WEAPON 

C  FWIL)  »  THE  FRACTION  BY  WEIGHT  OF  THE  ELEMENT  OCCUR I NG  IN  THE 
C  SOIL  OF  INTEREST 

C  HLII.L)  *  THE  HALF-LIFE  OF  DAUGHTER  ISOTOPE  OF  I  OF  L  IN  SECONDS* 

C  MINUTES*  HOURS*  DAYS*  OR  YEARS 

C  HOB  «  HEIGHT  OF  BURST  IFFETI 

C  ISO  •  ARRAY.  THE  NUMBER  OF  ISOTOPES  PER  ELEMENT  OCCUR ING  IN  NATURE 
C  KEV(I.L)  «  THE  NUMBER  OF  DIFFERENT  PHOTON  ENFRGIES  EMITTEO  BY  THE 
C  DAUGHTER  PRODUCED  BY  THE  ABSORPTION  OF  A  NEUTRON  BY  THE  ISOTOPE  I 

C  OF  THE  ELEMENT  L 

C  LH(I.L)  •  AN  INDEX  DENOTING  THE  UNITS  OF  THF  HALF-LIFE  OF  THE 
C  DAUGHTER  OF  I  ft-  .  Y ;  </  y,  S’*  tf*.) 

C  LMAC  *  NUMBER  OF  ELEMENTS  (UP  TO  20)  IN  CASING 

C  LMAX  ■  NUMBER  OF  ELEMENTS  (UP  TO  20)  IN  SOIL. 

C  NA  «  ARRAY.  THF  ATOMIC  MASS  NUMBER  OF  THE  ISOTOPE  I  OF  THE  ELEMENT  L 
C  NZ  •  ARRAY.  ATOMIC  NUMBER  OF  THE  ELEMENT 
C  RHOC  -  DENSITY  OF  CASING  IN  GM/CM**3 

C  RNY(N.I.L)  -  THF  EXPOSURE  RATE  WHICH  WOULD  BE  MEASURED  AT  3  FT  ABOVE 
C  AN  INFINITE  PLANE  UNIFORMLY  CONTAMINATED  WITH  A  SOURCE  EMITTING 

C  ONE  PHOTON  OF  A  PARTICULAR  ENERGY  PER  CM**2  PER  SEC. 

C  MR/HRI  /  PHOTON  /  (CM*«?-S€C>) 

C  SCHOB  *  SCALED  HEIGHT  OF  BURST 
C  SDRL(J)  *  R  /  HR  /  FISSION  /  CM**2  AT  TIME  J 

C  SIGC  »  THERMAL  ABSORPTION  CROSS  SECTION  OF  WEAPON  CASING  (BARNS) 

C  SIGH  I. L)  *  THF  THERMAL  NEUTRON  ABSORPTION  CROSS  SECTION  OF 

C  ISOTOPE  I  OF  L  (BARNS) 

C  SIGISC(I.L)  *  THE  THERMAL  NEUTRON  SCATTERING  CROSS  SECTION  OF  ISOTOPE 
C  I  OF  L  (BARNS) 

C  SIGH  «  MACROSCOPIC  ABSORPTION  CROSS  SECTION  OF  WEAPON  CASING  (CM**— 1 1 
C  SIGS  »  THERMAL  ABSORPTION  CROSS  SECTION  OF  SOIL  ( BARNS I 
C  SIGSSC  «  SCATTERING  CROSS  SECTION  OF  SOIL  FOR  THERMAL  NEUTRONS 
C  (BARNS) 

C  THICK  •  THICKNESS  OF  WEAPON  CASING  IN  CM. 

C  WYK  •  WEAPON  YIELD  IN  KT 
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cvj  oj  irj  vo|  t-j  *-j 


FORTRAN  INPUT  CARD  LIST  -  Program  Title  -  Induced  Activity  (p.  \  of  8) 
LIST  TITIZ 


CARD  WORD  COLS.  FORMAT  MNEMONIC  UNITS 

no.  of 

12  NTIMES  times 


1-10  F10.0 


01-20  F1Q.Q _ 1(2)  SSC, 

.21=30  nSLQ. _ IX3I-  .  sec. 

3LJifl  F1Q..Q _ 2(kl..  -jsec. 


fio.q _ 


51=60  £1£LH _ K61-  ...sea.. 

6l=2Q  F1Q.Q _ 1(7.)  -  sea. 


Number  of 

T mac  nn.nf  T.  casing  elements  (<  20] 


2? 


FORTRAN  INPUT  CARD  ARRAY  WORKSHEE 


FORTRAN  INPUT  CARD  ARRAY  WORKSHEET 


o 


FORTRAN  INPUT  CARD  ARRAY  WORKSHEET 


FORTRAN  INPUT  CARD  ARRAY  WORKSHEET 

♦Header  Card  Word  Colt.  Format  Mnemonic  Unite 

T  1^2”  “f?  LmaK  tfo.  of  toil  elementt  (LMAX  <  20) 

FORTRAN  INPUT  CARD  ARRAY  WORKSHEET  -  Program  Title- Induced  Activity  (P.5  of  8) 
ARRAY  TITLE  -  Induced  Activity  Soil  Element  Table 
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♦Foil ova  Induced  Activity  Casing  Isotope  Photon  Energy  Table  and 
Precedes  Induced  Activity  Soil  Element  Table. 
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FORTRAN  INPUT  CARD  ARRAY  WORKSHEET 


CARD  LIST 


FORTRAN  INPUT 


FORTRAN  INPUT  CARD  LIST  -  Progrm  Title  -Induced  Activity  (p.8  of  8) 
LIST  TITLE  BOMPAR  (Final  Card) 

CARD  WORD  COLS.  FORKAT  MNEMONIC  UNITS 

1  1-10  F10.0  WBC  JCT 


2  11-20  F10.( 


SOB 


ft. 


21-^0  Pio.o  THICK 
U-40  F10.0  RHOC 


Ul-50  F10.0  a^CEH  neutrons/flsslon 
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r\  n  n  n 


DIMENSION  T  ( 1 0  )  *  NZC20).  NZCC20)  •  1SOC70I.  FWC20I.  NAU0.20). 

2  AC10.70I.  KEVC10.201.  FAIC10.20I.  SIGIU0.20I.  HLC10.20)* 

i  LH| 10.201.  SAC  20 ) .  FAT  C  ?0 1 •  SFA I ( 20) .  FACC10.20).  XLAMC10.70). 

4  DRLU0.201,  FOGCS.l0.20j.  RNYCS.  10.201  •  FOGRNY 1 10.201 . 

5  OR  I  (  10.70!  .  (MCl  l'i.ZOI.  SORLHOI.  SORCCIOI.  SOKlOTi  lOi. 

6  0RKC10.70),  ORE C 10.20 J •  SDRKC10).  SDREC10). 

7  SIGISCC  10*20 )  .  SFAISCC20) 

KRO  *  5 
KTR  *  6 

read  input  tape  krd*  10. 

i  NTIMES 

READ  INPUT  TAPE  KRO.  12. 

1  CTIJ1.  J*l.  NTIMES1 

KK  *  1 
AC  •  0. 

STATEMENTS  100  TO  260  ARE  PFRFORMEO  TWICE.  THE  FIRST  TIME.  EXPOSURE 
RATE  FROM  ACTIVATED  ELEMENTS  OF  THE  BOMB  CASING  IS  CALCULATED. 

THE  SECOND  TIME.  EXPOSURE  RATE  DUF  TO  1N0UCE0  ACTIVITY  OF  THE  SOIL 
ELEMENTS  is  calculated. 

100  READ  INPUT  TAPE  KRD.  TO. 

1  lmax 

IF  IKK-11  101.  101.  102 

101  LMAC  *  LMAX 

102  READ  INPUT  TAPE  KRD.  14. 

1  CNZCL1 •  ISOCLI .  FWCLI.  L*l.  LMAX) 

DO  103  L  *  1.  LMAX 
IS  *  ISOCLI 

10)  READ  INPUT  TAPE  KRO.  16. 

1  IHACI.LI.  KEVII.LI.  FAIfl.LS.  SIGHT, L).  HLCI.LT.  LHCI.L). 

2  SIGISCC I.L) *  1*1*  IS) 

DO  104  L  *  1.  LMAX 
IS  *  ISOCLI 
DO  104  1*1.  IS 
KE  *  KEVCI.L) 

104  READ  INPUT  TAPE  KRO.  10. 

1  CFOGCN* I.L I .  RNYCN^I.Ll.  N=l.  KF1 

SIGS  =  0. 

SIGSSC  *  O. 

SFWMAX  =  0. 

DO  120  L  =  1 •  LMAX 
SA  CL  I  *  O. 

IS  =  ISOCLI 
DO  110  1=1.  IS 
AC  I .Li  =  NACI.LI 

110  SACL1  =  SACL1  *Ac I.L1«FAICI .LI 
120  SFWMAX  =  SFWMAX  ♦  FWCLl/SACLl 
DO  150  L  =  1,  LMAX 
FATIL1  =  FWCLI/CSACL»»5FWMAX1 
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sfakli  *  0. 

SFAJSC4LI  *  0 

IS  •  I SO* L* 

00  140  1  *  1*  IS 

SFAISC(L)  *SFA  I  SC  IL  !♦  FA1 d .Lt»SIGISCt 1 .LI 

140  SFAKU  »  SFAl(L)  ♦  FA! ( t .L t «SIGI ( 1 *L I 
IF  ncic-i»  141.  141*  1*9 

141  AC  -  AC  ♦  FATIU»SA<L> 

149  SIGSSC  *  SIGSSC  ♦  FAT(Lt*SFAISC(Lt 

150  S!GS  *  SIGS  ♦  FATILI*SFAItL) 

IF  (ICK-l)  151.  151.  152 

151  SIGC  *  SIGS 

152  00  200  L  *  1 •  LMAX 
IS  *  1  SOIL) 

00  200  I  «  1,  IS 

FACCI.L*  «  FATlLl»FAIII,L»*SIGItI.Ll/SIGS 
LG  *  LH(I.l) 

GO  TO  (200.140 .170 .180. 190 I •  LG 
160  Hl(I.L)  *  HL(1.L1*60. 

60  TO  200 

170  MLCI.L)  *  HLfl .LI *3600. 

GO  TO  200 

180  HL(l.L)  *  HL ( 1 «L 1 *86400. 

GO  TO  200 

190  HlC.LI  *  ML!  I  .L» •31556000. 

200  XLAMd.Lt  *  .695/HLII.L) 

00  201  J  *  1.  NTIMES 

201  SORL(J)  =  0.0 

00  260  J  *  1,  NTIMCS 
00  2*0  L  *  l.  LMAX 
ORLIJ.Li  =  O. 

IS  *  IS0IL1 
OO  220  T  *  1.  IS 
FOGRNYd.Lt  *  0. 

KE  =  KEVl I.Lt 
OO  210  M  *  1.  KZ 

210  FOGRNYd.Lt  =  FOGRNY(  I  »L1  ♦  FOGfN.I  »L)*RNY(N»I  «LI 

Will  1.1  1  *  FAC  1 1  *L  I  *XLAM(  I.Lt  •FOGRNY  (  I  *L  t*EXPF  (  **XLAH(  I«L)*TU)t 
220  DRL(J.L>  »  ORLU.Lt  ♦  0RKI.L1 
IF  (xiC-il  210.  230*  2*0 
230  ORClJ.Lt  *  ORLIJ.Lt 
N2CIL1  *  N7(l» 

240  SORLCJl  *  SORLCJt  ♦  ORLU.Lt 
IF  (KK-11  250.  250.  760 
250  SORCU!  *  SORLUI 
260  CONTINUE 

GO  TO  (270.280).  KK 

270  KK  =  2 

GO  to  100 

C  REDUCTION  OF  SOIL  EXPOSURE  RATE  IS  NOW  MADE  ACCORDING  TO  THE 

C  FRACTION  OF  NEUTRONS  ESCAPING  ABSORPTION  BY  THE  BOMB  CASING. 

C  THE  SOLID  ANGLE  FRACTION  OF  THE  FIREBALL  TOUCHING  THE  GROUND  AT 

C  THE  TIME  OF  HYDRODYNAMIC  SEPARATION.  AND  THE  FRACTION  OF 

C  NEUTRONS  THAT  REMAIN  »N  THE  SOIL 
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2*0  ALB  «  1.15*»«»SORTFISir,S/tSlGSSC  *  SIGS1I 

READ  INPUT  TAPE  KRO.  10. 

1  NWEP 

00  330  M  *  1  ,  NWEP 
READ  INPUT  TAPE  KRO*  12. 

1  WYK.  HOB*  THICK.  RHOC.  EMITM 


SCHOB  *  HOB/wYK**0.  13333333 

FOM  *  l.-SCHOB/SORTF(4.24*SCHOB*»2-234.*SCHOB«-4225.1 

SIGH  *  RH0C»0.A023*SIGC/AC 
EMI TC  *  EM!TN»CXPF|-S1GM*THICKI 
00  300  J  *  1.  NTIMES 
00  290  L  3  1 •  LMAX 
ORK(J.L)  *  ORCI J.L»*«f HITN-EMITCI 
290  DRE(J.L)  3  ORLt J.L»«EM!TC»FOM*ALB 
SDRKIJ)  3  SDRCI JMIEMITN-EMITC) 

SDREIJ*  *  SORLf J»»EM!TC*F0M«AL8 
300  SORTOT(J)  3  SORKIJ)  ♦  SORFiJ) 

WRITE  OUTPUT  TAPE  KTR.  30. 

1  WYK .  HOB.  SCHOB.  THICK.  AC.  RHOC*  EMITN.  EMITC.  SIGC.  SIGS. 

2  SIGSSC.  FOH.  ALB 

WRITE  OUTPUT  TAPE  KTR,  32. 

1  ITIJI.  J=1 ♦  NTIMES} 

00  310  L  3  i.  LMAC 
310  WRITE  OUTPUT  TAPE  KTR.  3*. 

1  NZCILI.  IORKI J.L* •  J*l.  NTIMES! 

WRITE  OUTPUT  TAPE  KTR.  36. 

1  ISORKtJI.  J31.  NTIMES* 

WRITE  OUTPUT  TAPE  KTR.  3S. 

1  ITIJI.  J31,  NTIMESI 

DO  320  L  3  1.  LMA* 

3?0  WRITE  OUTPUT  TAPE  KTR,  34, 

»  N2ILI.  IORFIJ.LI,  J=l,  NTIHESJ 

WRITE  OUTPUT  TAPE  KTR.  40, 

1  ISORFI JI.Jsl*  NTIMFSI 

330  WRITE  OUTPUT  TAPE  KTR,  42. 

1  ISORTOTIJI,  J*l.  NTIMFS1 


10  FORMAT 
12  FORMAT 
14  FORMAT 
16  FORMAT 
IB  FORMAT 


I I2.BX.F10. 01 
C7Fin,0l 
I12.IS.E13. 31 

II3.I2.F9.5.F1S.S.F10.3.I3.F15.3* 

IFIO.O.EIO.OI 
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-"BP"'"1- 


-  X 


30  FORMAT  I15M1  31/9  INDUCED  //  1?H  YfFLO  MCTl  •  F10.J.  / 

1  12H  HO«  IFEETI  «  FU«O,10X»11HSCALEO  HOft  .  F1J.3  ✓ 

2  7H  CASING*  /  6X,  ISHTHICKNrSS  JC.Mi  .F6.2.SX.12HAV.  AT.  WT.  • 

3  F7.3,  Sx*  20HAV.  DENSITY  IGM/CC)  »  F6*3*  / 

A  39H  N1IMRER  OF  NEUTRONS  EMITTED  PER  FISSION*  F3.2,  / 

5  A7H  NUMBER  OF  NEUTRONS  FSCApjNG  CASING  PFR  FISSION*  F5.2,  / 

6  31H  AVERAGE  CAPTURE  CROSS  SECTIONS*  /  6X*  9H0F  CASING. E10. 3. 

7  /  6X.7H0F  SOIL.E12.3/30H  SOIL  SCATTERING  CROSS  SECTION, E!0*3. 

8  /22H  SOLID  ANGLE  FRACTION  .FA. 4*  /  23H  SOIL  CAPTURE  FRACTION 

9  F6.4,  /✓/ I 

32  FORMAT  I10X,39HR/MR  PFR  FISSION/SOCM  DUE  TO  THE  CASING  // 

2  25X.  1 2H  TIME  ISFC)  / 

3  12H  ELFmFNT  121  •  5X.  I0F1O.P  /I 
3A  FORMAT  1 17,  lOX*  1P10E 10.31 

36  FORMAT  (12H0CASING  SUM  •  5X.  1PJ0E10.3  //» 

38  FORMAT  I ///SX.37HR/HR  PER  FISSION/SOCM  DUF  TO  THE  SOIL  // 

2  25X,  12H  TIME  ISfCl  / 

3  12H  ELEMENT  IZ1  *  SX.  10H0.0  /) 

AO  FORMAT  (12HO  SOIL  SUM  •  SX*  1P10E10.3  ) 

A2  FORMAT  (11K0CASING  AND*  / 

2  12H  SOIL  SUMMED*  SX.  1P10E10.3  ) 

CALL  EXIT 

END  (0*  0,  0.  0*  01 


38 


3179  INDUCE!) 


YIELD  IKT)  2.0 

MOB  « FEET)  -0. 

CALED  MOB 

-o. 

CASING 

THICKNESS  ICM)  -0. 

AV.  AT.  WT. 

36.415 

AV.  DENSITY  (GM/CC) 

NUMRF9  OF  NEUTRONS  EMITTED  PER  FISSION  1.00 
NUMBER  OF  NEUTRONS  ESCAPING  CASING  PER  FISSION  1.00 
AVERAGE  CAPTURE  CROSS  SECTIONS 
OF  CASING  0.103E+01 
OF  SOIL  0.203E-00 

SOIL  SCATTERING  CROSS  StCTION  0.916E+01 
SOLID  ANGLE  FRACTION  1.0000 
SOIL  CAPTURE  FRACTION  0.1700 


R/t-R  PER  FISS10N/S0CM  DUE  TO  THE  CASING 


TIME  ISECI 


ELEMENT  (2) 

3600.  4032. 

26 

0. 

0. 

13 

0. 

0. 

CASING  SUM 

0. 

0. 

R/HR  PER  FISSION/SOCM  DUE  TO  THE  SOIL 


TIME  I SEC) 


ELEMENT  (2) 

3600. 

4032- 

26 

3.742E-17 

3.742E-17 

64 

0. 

0- 

14 

9.019E-16 

8.738E-16 

19 

1.038E-14 

1.031E-14 

13 

2.529E-19 

2.889E-20 

11 

2.670F-13 

2.655E-13 

22 

2.247E-17 

9.505F-18 

1 

0. 

0. 

20 

2.196E-15 

1.221E-15 

17 

8.986E-15 

7.867E-15 

62 

0. 

0. 

25 

6.713E-13 

6.500E-13 

12 

2.256E-15 

1.334E-15 

5 

0. 

0. 

23 

0. 

0. 

24 

1.442E-18 

1.442E-18 

15 

0. 

0. 

8 

0. 

0. 

SOIL  SUM 

9.631E-13 

9.372E-13 

CASING  AND 
SOIL  SUMMED 


9.631E-13  9.372E-13 


APFEHDH  P 


Table  1  Elemental  Soil  Conpoaition 
Table  2  Isotope  Capture  parameters 
Table  3  Exposure  Bate  Paraaeters 
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TABIE  X 

Eleoent&l  Soil  Coapositicn 


Earth's 

Crust 

WES* 

(Jones) 

(Polan11)) 

Sleaent 

Z 

Weight 

Fraction 

(fA 

Weight 

Fraction 

(fw)k 

Weight 

Fraction 

H 

1 

1.4xlO"3 

-6 

0 

-4 

2.639BC10*2 

B 

5 

3.0x10 

1.0x10 

0 

0 

8 

4.66xlo"X 

5- 001x10* 1 

5. 47x10* 1 

Ka 

11 

2.83xl0“2 

i.34xio"2 

4.9xlO*3 

Mg 

12 

2.09xl0“2 

4. 4x10* 3 

2.09QO0"2 

Al 

13 

8.13xlo'2 

7 .63x10*  2 

4.8xlO”3 

Si 

14 

2-772x10" 1 

2.756x1c"1 

2.772X10*1 

P 

15 

I.l8xl0"3 

0 

l.lSxlO*3 

Cl 

17 

3.l4xlO_1* 

0 

3.l4xlO*4 

K 

19 

2.59X10"2 

2. 96x10* 2 

2.59Bd0"2 

Ca 

20 

'..63xl0“2 

8. 66x10* 2 

3.63x1c"2 

Ti 

22 

4.4xlO*3 

2.4xl0~3 

4.4xlO*3 

V 

23 

l.5xlo“4 

0 

0 

Cr 

24 

-4 

2.0x10 

0 

>4 

2.0x10 

Ha 

25 

l.QxlO"3 

3.0X10"11 

5-QxlO*11 

Fe 

26 

5- 0x10* 2 

1.12x10*  2 

5.0xl0"2 

So 

62 

6.5x10*^ 

6.5x10 

0 

Gd 

64 

6.4X10"6 

6.4x10*” 

0 

Modified  by  the  author  frcai  reference  2. 
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XABI£  2 


Isotope  Capture  Parameters'  ^  ^ 

Isotope 

fik 

Elemental 

Fraction 

i 

Absorption 
Cross  Section 
(barns) 

(*i)i 

Daughter 

Half-Life 

<°«>i 

Scat. 

Cross  Sect, 
(barns) 

H  1 

.99965 

.33 

1030** 

19 

H  2 

.00015 

.00057 

12.26y 

19 

B  10 

.1978 

3840 

1030. 

4 

B  11 

.8022 

.05 

.03s 

4 

0  1 6 

1.0000 

.0002 

1030 § 

4.2 

to  23 

1.0000 

.53 

15h 

4 

He  ak 

.7870 

.03 

1030 s 

3.6 

He  25 

.1013 

.27 

1030. 

3.6 

HS  a6 

.1117 

.03 

9-5n 

3.6 

Al  27 

1.0000 

.23 

2.3a 

1.4 

SI  26 

.9221 

.06 

1030 s 

1.7 

Si  29 

-0470 

.28 

1030 s 

1.7 

Si  30 

.0309 

.11 

2.62h 

1.7 

P  31 

1.0000 

.19 

l-3d 

5 

Cl  35 

.7553 

44.4 

3.08Kl05y 

16 

Cl  37 

.24*7 

.565 

37.5a 

16 

*  1030*  ■enu  the  daughter  is  stable. 


TABLE  2  (Cont.) 


Isotope 

fik 

Eleeental 

Fraction 

Absorption 
Cross  Section 
(bams) 

<Vi 

Daughter 

Half-Life 

Scat. 

Cross  Sect 
(bams) 

K  39 

•  9310 

2.2 

1.3*109y 

1.5 

K  40 

.0002 

70 

lO30* 

1.5 

K  4l 

.0688 

1.1 

12.4b 

1.5 

Ca  1*0 

.9 

.2 

1030 s 

3.2 

Ca  42 

.0064 

40 

lO30* 

3.2 

Ca  44 

.0206 

.7 

15*d 

3.2 

Ca  48 

.0018 

1.1 

8.5a 

3.2 

Ti  46 

-0793 

.6 

1030* 

4 

Ti  4? 

.0728 

1.7 

lO^s 

4 

Ti  48 

.7394 

8.0 

lO3^ s 

4 

Ti  4  ? 

.0551 

1.9 

1030* 

4 

Ti  50 

.0534 

.14 

5.8a 

4 

V  50 

.0024 

130 

1030 s 

5 

V  51 

.9976 

4.9 

3.76» 

5 

Cr  50 

.0431 

17.0 

27. 8d 

3 

Cr  52 

.8376 

.8 

I>. 

3 

Cr  53 

.0955 

18 

1030* 

3 

Cr  54 

.0238 

.38 

3.6a 

3 

Hn  55 

1.0000 

13.3 

2.58b 

2.3 
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TABX£  2  (Coot.) 


I 

Isotope 


Fraction 


Absorption 
Cross  Section 
(barns) 


<Vi  (°.c>l 

Daughter  Scat. 

Half-Life  Cross  Sect. 

(bams) 


Fs  5* 

.0582 

2.9 

2.96y 

11 

Fe  56 

.9166 

2.7 

1030* 

U 

Fe  57 

.0219 

2.5 

1030 8 

11 

Fe  58 

.0033 

1.1 

*6d 

11 

8a  1** 

.0316 

2 

*ood 

0 

Sm  1*7 

.1507 

87 

1030. 

0 

Sm  1*8 

.1127 

0 

1030 8 

0 

Sa  1*9 

.138* 

*0800 

1030 B 

0 

Sa  150 

.07*7 

0 

1030 8 

0 

Sa  152 

.2663 

22* 

*7h 

0 

sa  15* 

.2253 

5.5 

2*n 

0 

Gd  152 

.0020 

325 

230d 

0 

Gd  15* 

.0215 

0 

1030 a 

0 

Gd  155 

.1*73 

61000 

1030 8 

0 

Gd  156 

.20*7 

0 

1030 s 

0 

Gd  157 

.1568 

2*0000 

1030 8 

0 

Gd  158 

.2*187 

* 

l8.0h 

0 

Gd  160 

.2190 

.8 

3>6m 

0 

** 


TABLf*  3 

Exposure  Rate  Parameters 


Parent 

Isotope 

Ihoton 

Energy 

(Hey) 

*ViJ 

Photons  per 
Disintegration 

RJ 

photon/cn  -sec 
(xlO-6) 

M  1 

0 

0 

0 

H  2 

0 

0 

0 

B  10 

0 

0 

0 

B  11 

4.43 

.026 

15.13 

B  11 

3-22 

.013 

12.00 

0  1 6 

0 

0 

0 

Ha  23 

2.75 

1.000 

10.70 

Ha  23 

1-37 

1.000 

6.4o 

Ms  24 

0 

0 

c 

He  25 

0 

0 

0 

ME  26 

1.015 

.300 

5.05 

Mg  26 

.643 

.700 

4.33 

A1  2? 

1.78 

1.000 

7.75 

Si  26 

0 

0 

0 

Si  29 

0 

0 

0 

Si  30 

1.264 

.007 

6.04 

P  31 

0 

0 

0 

i 
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TABUS  3  (Coat.) 


Parent 

isotope 

Photon 

Energy 

(Mev 

Photons  per 
Disintegration 

BJ 

photon/ca 

Cl  35 

0 

0 

0 

Cl  37 

2.16 

shio 

8.92 

Cl  SI 

1.59 

.310 

7.14 

K  39 

1.46 

.110 

6.70 

K  <10 

0 

0 

0 

X  hi 

1.53 

.180 

6.95 

Ca  hC 

0 

0 

0 

Ca  42 

0 

0 

0 

Ca  44 

0 

0 

0 

Ca  48 

4.68 

.006 

15.76 

Ca  48 

4.C5 

.010 

14.10 

Ca  48 

3.10 

.890 

11.68 

Tt  46 

0 

0 

0 

T1  47 

Q 

0 

0 

T1  48 

0 

0 

0 

TL  49 

0 

0 

0 

Ti  50 

.928 

.042 

4.68 

Ti  50 

.605 

.013 

3.23 

Ti  50 

.323 

.945 

1.75 

V  50 

0 

0 

0 

V  51 

1.433 

1.000 

46 

6.62 

-»ec 


•UnT  T> 


\VWH»« 


\ 

# 


Parent 

Isotope 

Riaton 

Energy 

(Mev) 

Photons  per 
Disintegration 

.  HJ 

photon/cm 

(xKf6) 

Cr  50 

.323 

.090 

1.75 

Cr  52 

0 

0 

0 

Cr  53 

0 

0 

0 

Cr  54 

0 

0 

0 

Hn  55 

2.96 

.005 

11.34 

Mn  55 

2.65 

.018 

io.4o 

Mn  55 

2.13 

.195 

8.84 

Mr»  55 

1.81 

.262 

7.85 

Mi  55 

.845 

•977 

4.55 

Fe  54 

0 

0 

0 

Fe  56 

0 

0 

0 

Fe 

0 

0 

0 

Fe  58 

1.209 

.430 

6.10 

Fe  58 

1.191 

.030 

5.77 

Fe  58 

1.098 

.570 

5.38 

Sx  144 

0 

0 

0 

So  l4? 

0 

0 

0 

So  l48 

0 

0 

0 

So  149 

0 

0 

0 

3b  150 

0 

0 

0 

So  152 

0 

0 

0 

So  154 

0 

0 

4r 

0 

TABIE  3  (Coot.) 


Parent 

Isotope 

Photon 

Energy 

(Her) 

Rtotons  per 

Dia integration 

Jt/bt-- — % 

photon/ca  - 
(xlO~6) 

Gd  132 

0 

0 

0 

Gd  134 

0 

0 

0 

Gd  155 

0 

0 

0 

Gd  156 

0 

0 

0 

Gd  157 

0 

0 

0 

Gd  156 

0 

0 

0 

Gd  l60 

C 

0 

0 
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